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Abstract 

We investigate quantitatively the so-called leverage effect, which corre- 
sponds to a negative correlation between past returns and future volatility. 
For individual stocks, this correlation is moderate and decays exponentially 
over 50 days, while for stock indices, it is much stronger but decays faster. 
For individual stocks, the magnitude of this correlation has a universal 
value that can be rationalized in terms of a new 'retarded' model which 
interpolates between a purely additive and a purely multiplicative stochas- 
tic process. For stock indices a specific market panic phenomenon seems 
to be necessary to account for the observed amplitude of the effect. 



1 Introduction 

Several 'stylized fact' of financial markets, such as 'fat tails' in the distribution of 
returns or long ranged volatility correlations, have recently become the focus of 
detailed empirical study [p], |], ^ £|, [|, || [7|, §, |§ . Simple agent based models have 
been proposed, with some degree of success, to explain these features [[II], [II], |T2|, 



13, [L4L 19, |16fl. Another well-known stylized fact is the so-called 'leverage' effect, 



first discussed by Black |y| pL9| , who observed that the volatility of stocks tends 
to increase when the price drops. This effect is particularly important for option 
markets: not only does it imply that at-the-money volatilities tend to increase 
after price drops, but also that a significant skew in the volatility smile should 
appear |2l], [7], |(], |22| , as is indeed observed on markets where the leverage effect 
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is strong. This skew reflects the fact that a negative volatility-return correlation 



induces a negative skew in the distribution of price returns themselves [pT , pO 



Although widely discussed in the economic and econometric literature || |HJ , 
the leverage effect (or volatility-return correlation) has been less systematically 
investigated than the volatility clustering effect (volatility- volatility correlation). 
For example, one would like to know if the volatility-return correlation shows a 
long term dependence similar to that observed on the volatility-volatility corre- 
lation. Although various single correlation coefficients quantifying the leverage 
effect have been measured and discussed within garch models || [lTj, the full 
temporal structure of this correlation has not been quantitatively investigated. 
The economic interpretation of this leverage effect is furthermore still contro- 
versial; a recent survey of the different models can be found in W7\. Even the 



causality of the effect is debated fl7]] : is the volatility increase induced by the 
price drop or conversely do prices tend to drop after a volatility increase? 

In the present paper, we report some empirical study of this leverage effect 
both for individual stocks and for stock indices. We find unambiguously that 
correlations are between future volatilities and past price changes. For both stocks 
and stock indices, the volatility-return correlation is short ranged, with however 
a different decay time for stock indices (around 10 days) and for individual stocks 
(around 50 days). The amplitude of the correlation is also different, and is much 
stronger for indices than for stocks. We then argue that the leverage effect for 
stocks can be interpreted within a simple retarded model, where the absolute 
amplitude of price changes does not follow the price level instantaneously (as 
is assumed in most models of price changes, such as the geometric Brownian 
motion). Rather, absolute price changes are related to an average level of the past 
price. This reflects the lag with which market operators change their behavior 
(order volumes, bid-ask spreads, transaction costs, etc.) when the price evolves: 
the proportionality between absolute price changes and the level of the price 
itself is only ensured on the long run. We then show that this model does not 
represent adequately the leverage effect for stock indices, which seems to reflect 
a 'panic'-like effect, whereby large price drops of the market as a whole triggers 
a significant increase of activity. 



2 Empirical results 

We will call Si(t) the price of stock i at time t, and SSi(t) the (absolute) daily price 
change: 5Si(t) — Si(t + 1) — Si(t) . The relative price change will be denoted as 
Sxi(t) = 5Si(t)/ Si(t). The leverage correlation function which naturally appears 
in the calculation of the skewness of the distribution of price changes over a 
horizon T is [ 20 ]: 

C i (T) = U[Sx i (t + r)} 2 5x i (t)), (1) 
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Figure 1: Return- Volatility correlation for individual stocks. Data points are the 
empirical correlation averaged over 437 US stocks, the error bars are two sigma 
errors bars estimated from the inter-stock variability. The full line shows an 
exponential fit (Eq. (|3])) with A s = 1.9 and T s = 69 days. Note that £(r = 0) is 
not far from —2, as our retarded model predicts. 

which measures the correlation between price change at time t and a measure of 
the square volatility at time t + r. The coefficient Z is a normalization that we 
have chosen to be: 

Z = (fe,(t) 2 ) 2 (2) 

for reasons that will become clear below. (Note that since Sx is dimensionless, 
Ci(r) is also dimensionless, despite this apparent lack of homogeneity between 
the numerator and denominator.) 

In the following, we will consider r > 1. Negative values of r lead to very 
small values of the correlation function, indistinguishable from noise. In other 
words, the correlation exists between future volatilities and past price changes; 
conversely, volatility changes do not convey any useful information on future price 
changes. 

We have analyzed a set of 437 US stocks, constituent of the S&P 500 index 
and a set of 7 major international stock indices. Our dataset consisted of daily 
data ranging from Jan. 1990 to May 2000 for stocks and from Jan. 1990 to Oct. 
2000 for indices. We have computed Li both for individual stocks and stock 
indices. The raw results were rather noisy. We have therefore assumed that 
individual stocks behave similarly and averaged Li over the 437 different stocks 
in our dataset to give £5, and over 7 different indices (from the US, Europe and 



3 



10 



-10 



-20 



-30 




Empirical data 
Exponential fit 



50 



100 

T 



150 



200 



Figure 2: Return- Volatility correlation for stock indices. Data points are the 
empirical correlation averaged over 7 major stock indices, the error bars are two 
sigma errors bars estimated from the inter-index variability. The full line shows 
an exponential fit (Eq. with Aj = 18 and Tj = 9.3 days. 

Asia), to give £/. The results are given in Fig. 1 and Fig. 2 respectively. As 
can be seen from these figures, both C$ and £/ are clearly negative: price drops 
increase the volatility - this is the so-called leverage effect. These correlation 
functions can rather well be fitted by single exponentials: 



£s,/(r) = -A SJ exp 



T 



Tsj 



(3) 



For US stocks, we find As = 1.9 and T$ — 69 days, whereas for indices the 
amplitude Ai is significantly larger, Aj = 18 and the decay time shorter: Tj ~ 9 
days. This exponential decay should be contrasted with the very slow, power-law 
like decay of the volatility correlation function, which cannot be characterized by 
a unique decay time [0, || [5], || [7|, Therefore, a new time scale seems to be 

present in financial markets, intermediate between the very high frequency time 
scale seen on the correlation function of returns (several minutes) and the very 
low frequency time scales appearing in the volatility correlation function. 



3 A retarded model 

Traditional models of asset fluctuations postulate that price changes are propor- 
tional to prices themselves. The price increment is therefore written as: 



5Si(t) = S l (t)a i e(t), 



(4) 
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where Oi is the volatility and e a random variable with zero mean and unit 
variance, independent of all past history [^3] . Eq. (d) means that price increments 
are at any time proportional to the current value of the price. Although it is 
true that on the long run, price increments tend to be proportional to prices 
themselves, this is not reasonable on short time scales. Locally, prices evolve 
in discrete steps (ticks), following buy or sell orders that can only be expressed 
an integer number of contracts. The mechanisms leading to price changes are 
therefore not expected to vary continuously as prices evolve, but rather to adapt 
only progressively if prices are seen to rise (or decrease) significantly over a certain 
time window. The model we propose to describe this lagged response to price 
changes is to replace Si in Eq. @ by a moving average Sf 1 over a certain time 
window. More precisely, we write: 



6Si(t) = S«(t) a,e, S* (t) = £ £(r)Si(t - r), (5) 

T=0 

where /C(r) is a certain averaging kernel, normalized to one: 

oo 

£x:(t) = i. (6) 

r=0 

For example, an exponential moving average corresponds to K(t) = (1 — a)a T , 
(a < 1). It will be more congenial to rewrite S 1 / 2 as: 



Sf(*) = EK(r) 

r=0 



S l (t)-J2^(r')6S l (t-T'), (7) 



t'=1 



t'=1 

where /C(r) is the (discrete) integral of /C(r): 

oo 

X{t) = E /C(r'). (8) 

t'=t 

Note that from the normalization of /C(r), one has /C(0) = 1, independently of the 
specific shape of /C(r). This will turn out to be crucial in the following discussion. 

For the exponential model, one finds /C(r) = a T . The limit a — > 1 corresponds 
to the case where S^(t) is a constant, and therefore Eq. (]5|) corresponds to an 
additive model. The other limit a — ► (infinitely small averaging time window) 
leads to S^(t) = Si(t) and corresponds to a purely multiplicative model. An value 
of a close to one, a = 1 — e corresponds to an additive model on short time scales 
(<C T = e _1 ) and to a multiplicative model for long time scales (3> T) 0. A 
formulation of this model in terms of product of random matrices is given in the 
Appendix. 

In the following, we will assume that the relative difference between Si and 
Sf* is small. This is the case when: 



r) — a 



\ 



E ^V) « i (9) 
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For the exponential model, this is equivalent to ayT/2 < 1. In other words, 
i] ^ 1 means that the averaging takes place on a time scale such that the relative 
changes of price are not too large. Typically, T = 50 days, a = 2% per square- 
root day, so that 77 ~ 0.1. 

Let us calculate the correlation function Ci(r) to first order in 77 <C 1. One 
finds, using 5xi = 5 Si /Si and Eq. ([5]): 

[5 Xi (t + r)f Sxm) = <Z ( f 1 - 2 t nr') 6Sl ^ + T ~ T,) ) 6 4§-) ■ (10) 



S t (t + r) J Si(t) 

To first order in 77, one can replace in the above expression 5Si(t + r — t')/ Si(t + r) 
by (Tje(t + r — r') and SSi(t)/Si(t) by (Tie(t). Since the e are assumed to be 
independent, of zero mean and of unit variance, one has 



(e{t + T-T')e{t)) = 5 T!T , (11) 



Therefore: 



([^(t + r)] 2 Sxi(t)) = -2o\K{t) (12) 
With the chosen normalization for £j(t) (see Eq. (0)), we finally find: 

d(T) = -2X(t) (13) 

A very important prediction of this model is therefore that £j(r — * 0) = —2. As 
shown in Fig. 1, this is indeed rather well obeyed for individual stocks, with /C(r) 
given by a simple exponential. We have confirmed this finding by analyzing a 
set of 500 European stocks and 300 Japanese stocks, again in the period 1990- 
2000. The results are given in Figure 3 for the European stocks. We again find 
an exponential behavior with a time scale on the order of 40 days and, more 
importantly, and initial values of Ci close to the retarded model value —2 p4 |. A 



similar graph was obtained for Japanese stocks as well, which is interesting since 
this market did behave very differently both from the U.S. and European markets 
during the investigated time period. For the Japanese market, the prediction 
C(t — > 0) = — 2 not as good: an exponential fit of C(t) gives and As = 1.5 and 
Ts = 47 days. 

We therefore conclude that the leverage effect for stocks might not have a very 
deep economical significance, but can be assigned to a simple 'retarded' effect, 
where the change of prices are calibrated not on the instantaneous value of the 
price but on an exponential moving average of the price. On the other hand, 
as Fig. 2 reveals, the correlation function for indices has a much larger value for 
t = and the above interpretation cannot hold. We will discuss this in more 
details in the next section. 
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Figure 3: Return- Volatility correlation for European stocks. Data points are the 
empirical correlation averaged over 500 European stocks. The full line shows an 
exponential fit (Eq. @) with As = 1.96 and Ts = 38 days for European stocks. 
Note again that £ for small r is very close to the retarded model value —2. 

4 A one factor leverage model 

Figure 2 shows that (i) the 'leverage effect' for indices is much stronger than that 
appearing for individual stocks, but (ii) tends to decay to zero much faster with 
the lag t. This is at first sight puzzling, since the stock index is, by definition, 
an average over stocks. So one can wonder why the strong leverage effect for the 
index does not show up in stocks and conversely why the long time scale seen 
in stocks disappears in the index. We want to show in this section that these 
features can be rationalized by a simple one factor model, where the 'market 
factor' is responsible for the strong leverage effect. Let us therefore write the 
stock price increments as: 

6Si(t) = S*{t) \P4(t) + e^t)] , (14) 

where <f>(t) is the return factor common to all the stocks, (3i are some time in- 
dependent coefficients, and e« are the so-called idiosyncrasies, uncorrelated from 
stock to stock and from the common factor 0. The market index I(t) is defined 
as a certain weighted average of the stocks: 

N 

/(t) = y>^(*), (is) 



7 



where Wi are certain weights, of order 1/N. From linearity, one finds the same 
relation between the 'retarded' quantities: 



N 



I R (t)=Y, Wi S«(t). (16) 



i=l 



Neglecting terms of order 1/y/N, one finds, after summing Eq. ffiM over all i's: 

8I{t)=^I R {t)<t>{t), (17) 
where j3 is the weighted sum of the /3's: 



= ^ T;i , as) 



that we will take as a constant. One can always define (f)(t) such that this constant 
is unity, which is the choice we make in the following. Now, we postulate that 
there exists an index specific leverage effect, resulting from an increase of activity 
when the market as a whole goes down, reflecting a panic like effect. Therefore 
we write: 

(fit + r)^\ ~ (0 2 (t + r)0(t)) = -r(r) (19) 

In the following, we will indeed assume that V{r) is small, as the data suggests, 
and we neglect all mixed panic-retarded effects. To linear order in the correla- 
tions, one then finds the following index leverage effect: 



5I(t + r)] 2 5I(t)\ 



I(t) / 



-r(r)-2/C(r) U\t)) . (20) 



lit + r) 

Hence, the correctly normalized volatility-return correlation function reads: 

Cj(t) = -2X{t) - 7 (t) (21) 

with: 

r(r) = 7 (r)a / 4 (22) 

where a] = J ((f) 2 (t)) is the market volatility. Therefore, one explicitly sees that 
the slowly decaying part JC(t) should also appear in Cj{t). The amplitude of this 
retarded correlation (2/C(0) = 2) is however only 10% of the observed correlation 
(£j(0) = 18.). We have fitted the observed correlation for indices by a sum of 
two exponentials, with only the parameters of the 'fast' one left free, the slow one 
being fixed by fitting individual stocks. The resulting fit (not shown) was not 
significantly better (nor worse) than the single exponential fit of Fig. 2. Given 
the amount of noise in the data, it is difficult to prove or disprove the presence of 
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the slowly decaying correlation. Nevertheless, we argue that it should be present 
for reasons of consistency between the index and its constituents. 

Let us finally estimate the contribution of 7(7") to the individual stock leverage 
effect. A simple computation gives, to lowest order: 

£ 1 (t) = -21C(t)-[3?(-)\(t). (23) 



Since the market volatility 07 is a factor 3 smaller than the volatility of individual 



stocks Oi [25|, the prefactor in front of 7(7") is of order 1/100. Therefore, even 
if 7(7") is ten times larger than /C(r), the influence of the market leverage effect 
on individual stocks is effectively suppressed due to relatively large ratio between 
the stock volatility and the market volatility. Again, due to the amount of noise 
in the data, it is difficult to confirm directly the presence of the 7(7") contribution 
in Ci{r). However, since this contribution is small and decays relatively fast, we 
believe that its role for individual stocks can safely be neglected. 



5 Conclusion 

In this paper, we have investigated quantitatively the so-called leverage effect, 
which corresponds to a negative correlation between future volatility and past 
returns. We have found that this correlation is moderate and decays over a few 
months for individual stocks, and much stronger but decaying much faster for 
stock indices. In the case of individual stocks, we have found that the magni- 
tude of this correlation can be rationalized in terms of a retarded effect, which 
assumes that the reference price used to set the scale for price updates is not 
the instantaneous price but rather a moving average of the price over the last 
few months. This interpretation, supported by the data on U.S., European and 
Japanese stocks, appears to us rather likely, and defines an interesting class of 
stochastic processes, intermediate between purely additive (valid on short time 
scales) and purely multiplicative (relevant for long time scales), first advocated 
in 0. For stock indices, however, this interpretation breaks down and a specific 
market panic phenomenon seems to be responsible for the enhanced observed 
negative correlation between volatility and returns (and in turn to the strong 
skews observed on index option smiles). Interestingly, this effect appears to de- 
cay rather quickly, over a few days. As a simple one factor model shows, the 
two effects (retardation and panic) are not incompatible and should both be 
present in individual stocks and stock indices. However the relative amplitude of 
the retardation effect for indices and of the panic effect for stocks are both too 
small to be unambiguously detected with only stock or index data. Note that 
in both cases, the correlation between volatility and returns appears to decay 
exponentially, in strong contrast to volatility-volatility correlations which decay 
as a power-law. This power-law behavior has recently lead to the construction 
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of a beautiful multifractal model ||. It would be interesting, from a theoretical 
point of view, to generalize this model to account for a scale invariant leverage 
effect. Work in this direction is underway. 
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Appendix 

The exponential retarded model can be reformulated in terms of a product of 
random 2x2 matrices. Define the vector V(t) as the set S^t), S(t). One then 
has: 

V(t + 1) = M(t)V(t) (24) 

where M.(t) is a random matrix whose elements are Ain(t) = a, A4u{t) = 1 — a, 
•M-2i(t) = cre(t) and M.22(t) = 0. Using standard results on products of random 
matrices, one directly finds that the large time statistics of both components of 
V is log-normal. 
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